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CsJ , Abstract. A brief and incomplete review of known integrable and (quasi)-exactly-solvable 

quantum models with rational (meromorphic in Cartesian coordinates) potentials is given. 
All of them are characterized by (z) a discrete symmetry of the Hamiltonian, (ii) a number 
of polynomial eigenfunctions, (Hi) a factorization property for eigenfunctions, and admit 
(iv) the separation of the radial coordinate and, hence, the existence of the 2nd order integral, 
(w) an algebraic form in invariants of a discrete symmetry group (in space of orbits) . 
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1 Introduction 

(N 

>. 

■ In this paper we will make an attempt to overview our constructive knowledge about (quasi)- 
exactly-solvable potentials having a form of a meromorphic function in Cartesian coordinates. 
All these models have a discrete group of symmetry, admit separation of variable(s), possess an 
(in)finite set of polynomial eigenfunctions. They have an infinite discrete spectrum which is 

. linear in the quantum numbers. All of them are characterized by the presence of a hidden (Lie) 

algebraic structure. Each of them is a type of isospectral deformation of the isotropic harmonic 
oscillator. 

Let us consider the Hamiltonian = the Schrodinger operator 

- T— I ■ 

■ n = -A + V{x), xgR'*. 

The main problem of quantum mechanics is to solve the Schrodinger equation 

Ti'i'ix) = E'^ix), ^'(x) G L2(M'^), 

finding the spectrum (the energies E and eigenfunctions Since the Hamiltonian is an infinite- 
dimensional matrix, solving the Schrodinger equation is equivalent to diagonalizing the infinite- 
dimensional matrix. It is a transcendental problem: the characteristic polynomial is of infinite 
order and it has infinitely-many roots. Usually, we do not know how to make such a diagonal- 
ization exactly (explicitly) but we can ask: Do models exist for which the roots {energies), some 
of the them or all, can be found explicitly {exactly)? Such models do exist and we call them 
solvable. If all energies are known they are called exactly-solvable (ES), if only some number 
of them is known we call them quasi- exactly- solvable (QES). Surprisingly, all such models the 



*This paper is a contribution to the Special Issue "Symmetry, Separation, Super-integrability and Special 
Functions (S**)". The full collection is available at |http:/ /www.emis.de /journals/SlGM A/S4.html , 
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present author familiar with, are provided by integrable systems. The Hamiltonians of these 
models are of the form 

2 

in the exactly-solvable case and 

1 A ~2 2 W{n) + T f. , 4 
TiQES = -:tA + + ^4 + ar^ + 6r^ 

in the quasi-exactly-solvable case, where w, ojk, F are parameters, VK(r2) is a function on unit 
sphere and r is the radial coordinate. In both cases there exists the integral 

T=^c' + w{n), (1) 

where C is the angular momentum operator, due to the separation of variables in spherical 
coordinates. 

Now we consider some examples from the ones known so far. 

2 Solvable models 

2.1 Case 0{N) 

The Hamiltonian reads 

1 / 9^ ■? ')\ viv-l) 
^oiN) = 2 E (- 9^ + ^ ^« j + (2) 

or, in spherical coordinates, 
2r^ 9r V 

^=1/:^. (4) 

The Hamiltonian ([2]) is 0{N) symmetric. It describes a spherical-symmetric harmonic oscillator 
with a generalized centrifugal potential. Needless to say that the Hamiltonian T-Io{n) ^ 
commute. 

Thus, T has common eigenfunctions with the Hamiltonian T-Lo(n)- The spectrum can be imme- 
diately found explicitly, and all eigenfunctions are of the type 

where Y^iy{i}) is a A^-dimensional spherical harmonics, J-Y^iy{Q) = ■jY^£y(0,). The Hamil- 
tonian ([2]) describes an A^-dimensional harmonic oscillator with generalized centrifugal term. 

Substituting in ([3]) the operator T by its eigenvalue 7 and gauging away ^'o = r^e ~ we arrive 
at the Laguerre operator 

'•o(K) = iiorHnom - -Eo)'i'oU., = -2td^ + (2w-i-!^-e'j d,, (5) 



Prom Quantum Aj\i (Calogero) to -ff4 (Rational) Model 



3 



where Eq is the lowest energy and the parameter i is chosen in such a way as to remove singular 
term oc ^ in the potential in ([3]). is the algebraic form of the Hamiltonian ([3]). The 
gauge-rotated Hamiltonian /lo(Ar) © is s/(2)-Lie-algebraic (see below), it has infinitely-many 
finite-dimensional invariant subspaces in polynomials Vn, n = 0, 1, . . . forming the infinite flag 
(see below), its eigenfunctions Pni'f'^ = i) are nothing but the associated Laguerre polynomials. 
By adding to ho(N) dS]) the operator 

J/j(qes) ^ 4 ^^^2 - 7) ^ - 4akt + 2ujk, (6) 

we get the operator ho(N) + Sh^'^'^'^^ which has a single finite-dimensional invariant subspace 
Vk = {t''\0<p<k), 

of the dimension (k + 1). Hence, this operator is quasi-exactly-solvable. Making the change of 

variable t = r'^ and gauge rotation with ^0 = t'^e~^ ~ we arrive at the 0(A^)-symmetric QES 
Hamiltonian [2] 



where 7, F, oj are parameters and 7 is replaced by the operator J^. In ^ a finite number of the 
eigenfunctions is of the form 

Pfe(r2)r2%^(J])e-'^-^, 

they can be found algebraically. It is worth noting that at a = the operator ho(N) + Sh^'^'^^^ 
remains exactly-solvable, it preserves the infinite flag of polynomials V and the emerging Hamil- 
tonian has a form of ([3]). 

2.2 Case (Za)^ 
The Hamiltonian reads 

t=i ^ ' ^ 1=1 

or, in spherical coordinates. 



where 



1 / \ 



i=l 

and J-" is given by ([5]). The Hamiltonian ([8]) is {1'2)^ symmetric. It defines the so-called Smoro- 
dinsky-Winternitz integrable system [1] which is in reality the maximally-superintegrable (there 
exist (2A^ — 1) integrals including the Hamiltonian) and exactly-solvable. Gauging away in ([8|) 

the ground state, 'I'o = 11 i^i)^ ^xp ( — 2^), and changing variables to tj = xf we arrive at the 

i=l 

algebraic form. Also it admits QES extension. The system described by the Hamiltonian ([S]) at 
Vi = v is a particular case of the i?CAr-rational system (see below). 
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<^ X2 <C X^ <^ — <C -^yv 

Figure 1. iV-body Calogero model. 

2.3 Case Ajv-i 

This is the celebrated Calogero model (Atv-i rational model) which was found in [3]. It descri- 
bes N identical particles on a line (see Fig. [1]) with singular pairwise interaction. 
The Hamiltonian is 

where the singular part of the potential can be written as 

2^ (^._r.)2 = 72 ' Wa^-i(^) = 2^[-, 3-1 , 

i>i ^ i=l Vt - T/ 



(10) 



Here r is the radial coordinate in the space of relative coordinates (see below for a definition) 
and WAff_^i^) is a function on the unit sphere. 

Symmetry: Sn (permutations Xi — )• xj) plus Z2 (all Xi — )■ —Xi). The ground state of the 
Hamiltonian ^ reads 

^o{x) = '[[\x,-Xj\''e~^^''h (11) 
Let us make the gauge rotation 

and introduce center-of-mass variables 

Y = ^Xi, yi = Xi-^Y, i = l,...,N, 
and then permutationally-symmetric, translationally-invariant variables 

{xi,X2, . . . ,xn) {Y,tnix) = an{y{x)) | ?i = 2, 3, . . . , 7V) , 

where 

(^k{x)= ^ Xi-^Xi^- ■ ■ Xi^, ak{-x) = {-)^ak{x), 

ii<i2<---<ik 

are elementary symmetric polynomials, and 



ii = 0, t2 y^jxj - Xj) 



2 2 
— r , 



i<j 



hence, the variable t2, which plays fundamental role, is defined by radius in space of relative 
coordinates. After the center-of-mass separation, the gauge rotated Hamiltonian takes the alge- 
braic form [4] 
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where 

(iV-i + l)(l-j) \- rn, 

■Aij = — ti^itj^i + 2^ [21 - J + 

l>max{l,j~i) 

Bi = ^(1 + vN){N - i + 2){N - i + l)ti_2 + 2uj{i - 1)U. 
Eigenvalues of ()12p are 

TV 
1=2 

hence, the spectrum is linear in the quantum numbers P2,3,...,N = 0,1,..., it corresponds to 
anisotropic harmonic oscillator with frequency ratios 1 : 2 : 3 : • • • : (A^ — 1). 

It is easy to check that the gauge-rotated Hamiltonian /icai has infinitely many finite- 
dimensional invariant subspaces 

= {t.PHsP-' ■ ■ ■ tiv^^ I < < n). 

where n = 0, 1, 2, . . .. As a function of n the spaces Vn^ form the infinite flag (see below). 



2.3.1 The grZrf+i-algebra acting by 1st order differential operators in M'* 

It can be checked by the direct calculation that the gld+i algebra realized by the first order 
differential operators acting in in the representation given by the Young tableaux as a row 
(n, 0, 0, . . . , 0) has a form 

d-l 

yJi ~ ^ = 1, 2 . . . , d, iTij = , z, J = 1, 2, . . . , d, = ^ ^ 

=tiJ^ =tii^tj^-T^, i = l,2...,d. (13) 

where n is an arbitrary number. The total number of generators is (d + 1)^. If n takes the 
integer values, n = 0, 1, 2, . . ., the finite-dimensional irreps occur 

It is a common invariant subspace for ([13]). The spaces P„ at n = 0, 1, 2, . . . can be ordered 



rocricr2c ■■■ cVnC ■■■ cv. 



Such a nested construction is called infinite flag {filtration) V. It is worth noting that the 
flag 'P^'^-' is made out of finite-dimensional irreducible representation spaces vif'^ of the algebra 
gld+i taken in realization (jl3p . It is evident that any operator made out of generators (jl3p 
has finite- dimensional invariant subspace which is finite- dimensional irreducible representation 
space. 
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2.3.2 Algebraic properties of the Calogero model 



It seems evident that the Hamiltonian (jl2p has to have a representation as a second order 
polynomial in generators (jl3p at d = N — 1 acting in M^~^, 

where the raising generators J'^^ are absent. Thus, gl{N) (or, strictly speaking, its maximal 
affine subalgebra) is the hidden algebra of the A^-body Calogero model. Hence, /icai is an 
element of the universal enveloping algebra UgK^j^iy The eigenfunctions of the A-body Calogero 

model are elements of the flag of polynomials 'p(^~^\ Each subspace is represented 

by the Newton polytope (pyramid). It contains C^_^^__-^ eigenfunctions, which is equal to the 
volume of the Newton polytope. 

Making the gauge rotation of the integral ([1]) with WAff_i{^) given by (jlOp 



/Cal = ^0 ^-^Cal - i^0)^0, 

where Fq is the lowest eigenvalue of the integral, /"cai^o = -^o^Oi the integral gets the algebraic 
form 

92 d 

where fij is 2nd degree polynomial in t, f2j = 0, and gi is 1st degree polynomial in t, g2 = 0. It 
also can be rewritten as the second degree polynomial in the gl{N) generators, 

/Cal = Pol2(Ji",j;/). 

2.3.3 sZ(2)-quasi-exactly-solvable generalization of the Calogero model 

By adding to /icai (fT2]) . the operator 



^/,(qes) ^ 4/^2 _ 4^^^2 + 2cok, 

^ ' dt2 

we get the operator hcai + Sh^"^^^"^ having finite-dimensional invariant subspace 
Vk = (tf |0 <p< A;). 

By making a gauge rotation of ^cai + dh^'^'^^^ and changing of variables to Cartesian ones we 
arrive at the Hamiltonian [S] 

1=1 ^ I -' j^i J' 

+ a^r'^ + 2aujr^ - a[2k + 2n{l + v + vn) - 7 - 1] r^. 

For the Hamiltonian, {k + 1) eigenfunctions are of the form 



n 



^(q-)(x) = n - XjV{ryPk{r^) exp 



„2 "^4 



2^ ' 4 

i=l 



(r2)^P,(r2)exp(-^r4) ^0, 
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where ^'o is given by ([TT]) . Pk is a polynomial of degree A; in = ^( r = t2. All 

i<j 

remaining eigenfunctions can be represented in the same form but P^'s are not polynomials 
anymore being functions depending on all variables Xj. It is worth noting that at a = the 
operator /icai + Sh^'^'^^^ remains exactly-solvable, it preserves the flag of polynomials 'p(^~^) and 
the emerging Hamiltonian has a form of ([9]) with the extra term ^ in the potential. Its ground 
state eigenfunction is (r^)'^^'o. It is the exactly-solvable generalization of the Calogero model ([9]) 
with the Weyl group W(A^_-^) as the discrete symmetry group, 

r 

2.4 Case: Hamiltonian reduction method 

In this methocj^ a family of integrable and exactly-solvable Hamiltonians associated with Weyl 
(Coxeter) symmetry was found with the Calogero model as one of its representatives. The idea 
of the method is beautiful and sufficient transparent, 

• Take a simple group G, 

• Define the Laplace-Beltrami (invariant) operator on its symmetric space (free/harmonic 
oscillator motion), 

• Radial part of Laplace-Beltrami operator is the Olshanetsky-Perelomov Hamiltonian re- 
levant from physical point of view. The emerging Hamiltonian is the Weyl-symmetric, it 
can be associated with root system, it is integrable with integrals given by the invariant 
operators of higher than two orders with a property of solvability. 

Rational case. This case appears when the coordinates of the symmetric space are intro- 
duced in such a way that the zero-curvature surface occurs. Emerging the Calogero-Moser- 
Sutherland-Olshanetsky-Perelomov Hamiltonian in the Cartesian coordinates has the form. 



1 ^ 
k=l 



U 2 2 



where i?+ is a set of positive roots, x is a position vector and i/j^i are coupling constants 
(parameters) which depend on the root length. If roots are of the same length, then v^^^ have 
to be equal, if all roots are of the same length like for A„, then all = v. In the Hamiltonian 
Reduction the parameters u\a\ take a set of discrete values, however, they can be generalized 
to any real value without loosing a property of integrability as well as of solvability with the 
only constraint of the existence of L^-solutions of the corresponding Schrodinger equation. The 
configuration space for (jl4p is the Weyl chamber. The ground state wave function is written 
explicitly, 

*o(^)= n l(«-^)r'"'e-""'/'- (15) 

The Hamiltonian p4p is completely-integrable: there exists a commutative algebra of inte- 
grals (including the Hamiltonian) of dimension which is equal to the dimension of the configura- 
tion space (for integrals, see Oshima [7] with explicit forms of those). For each Hamiltonian (jl4p 
after separation of center-of-mass coordinate (if applicable) the radial coordinate (in the space 
of relative coordinates) can be also separated. It gives rise to the existence of one more integral 



^For review and references see e.g. [6]. 
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of the second order ([T]). Hence, the Hamiltonian (|14p is super-integrable. The Hamiltonian ([T] 
is invariant with respect to the Weyl (Coxeter) group transformation, which is the discrete 
symmetry group of the corresponding root space. 

The Hamiltonian (|14p has a hidden (Lie)-algebraic structure. In order to reveal it we need 

to 

• Gauge away the ground state eigenfunction making similarity transformation {'^q)~^{'H — 
-Eo)^o = h, 

• Consider the Hamiltonian in the space of orbits of Weyl (Coxeter) group by taking the 
Weyl (Coxeter) polynomial invariants as new coordinates, these invariants are 

where o's are the degrees of the Weyl (Coxeter) group, 0, is an orbit. 

The invariants t are defined ambiguously, up to invariants of lower degrees, they depend on 
chosen orbit. 



2.5 Case BCn 

The i?CAr-rational model is defined by the Hamiltonian, 

N / ^2 



1=1 ' 



2 2 



i<j 



1 1 

+ 



1^2(1^2 - 1) J_ 



where u, u, V2 are parameters. If = 0, the Hamiltonian (jl6p is reduced to ([8] 
of the system is S'at ® (^2)^ (permutations xi — ?• Xj and Xi — ?• — Xj). 
The ground state function for ()16p reads 



(16) 

I. The symmetry 



N 



n 



11^2 



(cf. (lisp ). Making the gauge rotation 

hBC^ = (^0)"^ i^BC^ - Eo)^o, 
and changing variables 

{xi,X2, . . . ,xn) (o-fe(x2) I /c = 1,2, . . . ,iV) , 

where 



2 2 



ii<i2<---<ik 

where r is radius, we arrive at [S 

52 



cri(x^) = X^ + x| + • • 



I 2 _ 2 

• + Xj^ — r , 



+ S^(<^)7^, 
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with coefficients 

Aij = -2^(2/ + 1 + i - i)ai-i-iaj+i, 

l>0 

Bi = [1 + 1/2 + 2v{N -i)]{N -i + l)a,-i + 2ujiai. 

This is the algebraic form of the BCn Hamiltonian. Assuming polynomiality of the eigenfunc- 
tions we find the eigenvalues: 



N 
1 



hence, the spectrum is equidistant, linear in the quantum numbers and corresponds to anisotropic 
harmonic oscillator with frequency ratios 1 : 2 : 3 : • • • : A^. The Hamiltonian hscN infinitely 
many finite-dimensional invariant subspaces of the form 

= {aiP^a2^' ■ ■ ■ aN^^ \0<^p,<n), 

where n = 0, 1, 2, . . .. They naturally form the flag 'p(^). The Hamiltonian can be immediately 
rewritten in terms of generators as a polynomial of the second degree, 

hBC^=Poh{J.r,J^J^), 

where the raising generators J'^ are absent. Hence, gl{N + 1) is the hidden algebra of the -BC^r 
rational model, the same algebra as for the Ajv-rational model. The eigenfunctions of the BC^- 
rational model are elements of the flag of polynomials p(^). Each subspace contains C^_^_^ 

eigenfunctions (volume of the Newton polytope (pyramid) Pn^^). 

The BCjsf Hamiltonian admits 2nd order integral as result of separation of radial variable 

Evidently, the commutator 
Gauge-rotated integral 

where J^bCn'^o — ^o^o, takes the algebraic form in t-coordinates, 

where fij is 2nd degree polynomial, fij = 0, and Qi is 1st degree polynomial, gi = 0, 

fBC^=Poh{j.r,Ji,'>), 

in terms of the gl{N + 1) generators. It is worth mentioning that the commutator of [h,f] 
vanishes only in the realization (jl3|) . otherwise, 

[hBC^{J)jBC^{J)]^0. 
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2.5.1 si(2)-quasi-exactly-solvable generalization of the BCn rational model 

By adding to ^bcjv the operator 

J/,{qcs) ^ 4/ 2 _ 4q^^^ ^ 2i0k, 

which is the similar to one for the Calogero model, we get the operator HbCn + ^h^""^^^ which 
has the finite-dimensional invariant subspace 

Making a gauge rotation of Iibcn + 5h^^^^^ and changing the variables it's back to the Cartesian 
ones the Hamiltonian becomes 



1=1 ^ « ^ i<j 



+ 



Mi^2 - 1) A 1 27 [7 - 2jV(l + 2u{N - 1) + i^a) + 3] 
+ 2 ^ x2 + r2 

+ aV^ + 2awr'^ - a [2A; + 2iV(l + 2u{N - 1) + 2^2) - 7 - 1] r^ 

for which (k + 1) eigenfunctions are of the form 



n 



=X{\xi- x]\'J{\xr{ryPk{r^) exp 

i<j 1=1 

where is a polynomial of degree A; in = xj. 

i=l 

It is worth noting that at a = the operator Kbcn + 5h^'^^^^ remains exactly-solvable, it 
preserves the flag of polynomials V^^^ and the emerging Hamiltonian has a form of (jl6p with 
the extra term ^ in the potential. Its ground state eigenfunction is {r'^)'^'^^. It is the exactly- 
solvable generalization of the i?CAr-rational model (jl6p with the Weyl group W(BCjy) as the 
discrete symmetry group, 

r 

^W{BC^) =^BCjv + ^• 

Now we are in a position to draw an intermediate conclusion about and BCn rational 
models. 

• Both An- and -BCAr-rational (and trigonometric) models possess algebraic forms associated 
with preservation of the same flag of polynomials 'p(^). The flag is invariant w.r.t. linear 
transformations in space of orbits 1 1— )• t+A. It preserves the algebraic form of Hamiltonian. 

• Their Hamiltonians (as well as higher integrals) can be written in the Lie-algebraic form 

h = Voh{j{b(lglN+i)), 

where P0I2 is a polynomial of 2nd degree in generators J of the maximal affine subalgebra 
of the algebra h of the algebra qIn+i in realization ()13p . Hence, qIn+i is their hidden 
algebra. From this viewpoint all four models are different faces of a single model. 

• Super symmetric A^- and BCn -rational [and trigonometric) models possess algebraic 
forms, preserve the same flag of (super) polynomials and their hidden algebra is the super- 
algebra gl{N + 1|A^) (see |8]). 
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In a connection to flags of polynomials we introduce a notion 'characteristic vector'. Let us 
consider a flag made out of "triangular" linear space of polynomials 



In lattice space V 



= (xf ■■■xl'\^< fm + /2P2 + • • • + fdPd < n) , 

where the "grades" /'s are positive integer numbers and n = 0,1,2, 
defines a Newton pyramid. 

Definition 1. Characteristic vector is a vector with components /j: 
/ = (/l) /2, • • • , fd)- 

From geometrical point of view / is normal vector to the base of the Newton pyramid. The 
characteristic vector for flag "P^*^) is 

/o = (l,l,...,l) 



2.6 Case G2 
Take the Hamiltonian 

3 . 



i=l ^ ' i<j 



+ 3/x(//-l) 



{Xk +Xl- 2Xm)'^' 



(17) 



where uj, v, /i are parameters. It describes the Wolfes model of three-body interacting system [9] 
or, in the Hamiltonian reduction nomenclature, the G2-rational model. The symmetry of the 
model is dihedral group Dq. The ground state function is 



n 



3 1 V— * 2 

-i-r -9^ z2 



i<j k<l,k,ly^m 

Making the gauge rotation 

he, = {^orHnc, - E)^o, 
and changing variables 

Y = Yxi, yi = Xi-^Y, i = 1,2,3, {xi,X2,X3) 



where 



-Vi - 2/2 - ym ~ -r^, A2 = [yiy2{yi + y2)f, 



and separating the center-of-mass coordinate we arrive at 



+ {2ujXi + 2[l + 3{fi + u)]}dx, + 



6oj\2 - -{l + 2fi)\l 



9x2, 
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which is the algebraic form of the Wolfes model. The eigenvalues of Hq^ are 
e{p} = 2a; (Pi +3p2)- 

It coincides to the spectrum of anisotropic harmonic oscillator with frequency ratio 1 : 3. 

Separating the center-of-mass in (|17p and introducing the polar coordinates {g, ip) in the 
space of relative coordinates we arrive at the Hamiltonian 

/ N o9 1 o 1 o9 9 9 9v(u — 1) 9ii(u — 1) 

r ^ cos^ i^p sm 3(/3 

It is evident that the integral of motion which appears due to separation of variables in polar 
coordinates (cf. ([T|)) has the form 



^ cos^ sm^ 3(/9 



It is evident that after gauge rotation with ^-iid change of variables to (Ai, A2) the integral J- 
takes algebraic form. 

The Hamiltonian /iqj has infinitely many finite-dimensional invariant subspaces 

^Si,2) = (^f^? I < Pi + 2p2 < n), n = 0, 1, 2, . . . , 

hence the flag ^(12) with the characteristic vector / = (1,2) is preserved by h^^. The eigen- 

(2) (2) 
functions of /igj are are elements of the flag of polynomials V)^-^ 2) • Each subspace ^-^ 2) ~ 

(2) 

^n-i (1 2) contains ~n eigenfunctions which is equal to length of the Newton line Cnn = 

n{XiP'X2P''\pi + 2p2 = n). 

A natural question to ask: What about hidden algebra? Namely: Does algebra exist for which 
(1 2) ^'^ space of (irreducible) representation? Surprisingly, this algebra exists and it is, in 

fact, known. 

Let us consider the Lie algebra spanned by seven generators 
J^ = du Jn = tdt-'^, 4 = 2u9„-^, 

J^ = t^dt + 2tudu-nt, R^ = fdu, i = 0,l,2, L = {Rq, Ri, R2). (19) 

It is non-semi-simple algebra gl{2,'K) x 7^^^^ (S. Lie [lOj at n = and A. Gonzalez-Lopez et 
al. [llj at n 7^ (case 24)). If the parameter n in ()19p is a non-negative integer, it has 

= (tPu^jo < (p + 2g) < n), 
as common (reducible) invariant subspace. By adding 
= udl 

to gl{2,R) X 7^(2) (see (fUD), the action on V^^j^^ 2) 

gets irreducible. Multiple commutators of 
with Tq^'' generate new operators acting on V^J-^ 2)1 

rf ) ^ [j\ [j\ [. . . j\ rf ] . . .] = n92-Vo( Jo + 1) • • • ( Jo + i - 1), i = 0, 1, 2, 
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where Jq = tdt + 2udu — n, and all of them are of degree 2. These new generators have a property 
of nilpotency, 

rf ) = 0, i>2, 

and commutativity: 

[rf ) , rf = 0, i,j = 0,1,2, u= {4^^ , ) , rf . (20) 

plus (pUj) span a linear space with a property of decomposition: gi^^) = L x ((7/2 ® Jo) 
(see Fig. [2]). 

L , t , U 

Figure 2. Triangular diagram relating the subalgebras L, U and g£2- ^2(5^2) is a polynomial of the 
2nd degree in 17^2 generators. It is a generalization of the Gauss decomposition for semi-simple algebras. 

Eventually, infinite- dimensional, eleven- generated algebra {by (jl9p and Jq plus (j20p . so that 
the eight generators are the 1st order and three generators are of the 2nd order differential 
operators) occurs. The Hamiltonian /igj can be rewritten in terms of the generators ()19p . (j20p 
with the absence of the highest weight generator J^, 

hG2 = {J^ + ^J^)J^ - \j'^R2 + 2[3(/i + + 1] + 2wj2 + 3a; J3 - ^(1 + 2^j)R2, 

where J^'^ = Jq'^- Hence, g/(2,M) x TZ^"^^ is the hidden algebra of the Wolfes model. 

{i) G2 Hamiltonian admits two mutually- non-commuting integrals: of 2nd order as the result 
of the separation of radial variable (see (|18p ) and of the 6th order. If a; = the latter 
integral degenerates to the 3rd order integral (the square root can be calculated in closed 
form) . 

[ii) Both integrals after gauge rotation with ^0 take in variables Ai^2 the algebraic form. Both 
preserve the same flag "P^f 2) ■ 

{Hi) Both integrals can be rewritten in term of generators of the algebra g^'^^ : integral of 2nd 
order in terms of gl{2,M) ix TiP'^ generators only and while one of the 6th order contains 
generators from £ as well jl3j . 

2.6.1 sZ(2)-quasi-exactly-solvable generalization 

By adding to /1G2) the operator (the same as for the Calogero and the BCn models) 

0X1 

we get the operator hcz + Sh^"^^^^ having single finite-dimensional invariant subspace 
Vk = {Xl\0<p<k). 

Making a gauge rotation of /igj + Sh^'^'^^\ changing of variables (y, Ai^2) back to the Cartesian 
coordinates and adding the center-of-mass the Hamiltonian becomes 

1=1 ' ' i<j ■' 
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+ 3M/i-l) Yl 



+ 



47(7 + 3fi + 3u) 



+ a'^r^ + 2aur^ + 2a [2k - 3(^ + z^) - 2(7 + 1)] r^, 
for which {k + 1) eigenfunctions are of the form 



(qes) 



i<j 



n 



|xi + Xj - 2xp\^{ryPk{r^) exp 



'2 



'i=i 



where Pk is a polynomial of degree /c in r^. 

It is worth noting that at a = the operator /i^j + Sh^'^'^^^ remains exactly-solvable, it 
preserves the flag of polynomials ^(12) emerging Hamiltonian has a form of (jl7p with 

the extra term ^ in the potential. Its ground state eigenfunction is (r^)'''\I'o. It is the exactly- 
solvable generalization of the G2-rational model (jl7p with the Weyl group W(G2) as the discrete 
symmetry group, 



W(G2) 



+ 4- 



2.7 Cases F4 and -Eg, 7,8 

In some details these four cases are described in [12]. 

2.8 Case l2{k) 



In some details this case is described in [13]. It is worth noting that although the Hamiltonian 
Reduction nomenclature is assigned to this case the parameter k takes any real value. Discrete 
symmetry group of the Hamiltonian appears for integer k. 



2.9 Case ifg 

The H2, rational Hamiltonian reads 



1 ^ 



fe=i 



d"^ no Uiv - 1) 



+ 2V(V-1) V V ^ ^ rrr, (21) 

where {i, j, A;} = {1, 2, 3} and all even permutations, w, v are parameters and 
1± ^5 

the golden ratio and its algebraic conjugate. Symmetry of the Hamiltonian (|2ip is the H3 
Coxeter group (the full symmetry group of the icosahedron) . It has the order 120. In total, the 
Hamiltonian (j2ip is symmetric with respect to the transformation 



Xi < > Xj, (p+ i > (f^. 

The ground state is given by 







A'(A^2 exp 
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where 

3 

Ai = j]xfc, n n [xi + i-ir^+xj + i-ir^if.xk]. 

k=l {jj'.fe} A*!, 2=0,1 

Making the gauge rotation 

we arrive at new spectral problem 

hH-j,(t>{x) = -2e(j){x). 

After changing variables (3:1^2,3 ''"1,2,3): 

2,2,2 2 
Ti = Xi+ X2 + = r , 

r2 = -^{xl + xl + xl) + ^(2 - 5ip+){xlxi + xlxi + xlxf) 

3 . \/24 24 2 4\ 

+ Y^(2 - 599_)(a;iX3 + x^x^ + XgXa) - —, 
T3 = + xl^ + 40) + ^(1 + 5(^-)(xfxi + xlxl + x«x?) 

+ ^(1 + 5(^+)(x?xi + x«x? + xlxl) + ^(1 - 5<^-)(x?x^ + xlxi + x^xf) 

+ 25 ^''^ ~ ^9^+) (^1^3 + ^2^1 + ^3^2) ^(^1^2^3 + X%x\x\ + X3X^X2) 

( ry* ^ ry* ^ ry*^ I ry* ^ ry*^ ^Y*^ I ry* ^ ry*^ 

25 l,-'^l'^2-^3 "T ■''2-'^3''^l "T •^3-^l-^2j ) 

in the gauge-rotated Hamiltonian, it emerges in the algebraic form [14j 

t,]=i -I j=i 

where 

All = 4ri, A12 = 12r2, ^13 = 2OT3, 

48 2 45 . 16 2 2 . 64 128 o 

-422 = -Ti T2 + — r3, A23 = TTnT2 - 24ri r3, ^33 = -T1T2TS + —T2 , 

5 2 15 o 45 

48 

= 6 + 60Z^ - 4WT1, ^2 = (1 + 5iy)T? - 12UJT2, 

5 

64 

53 = - — (2 + 5i/)Tir2-20wr3, 
15 

which is amazingly simple comparing the quite complicated and lengthy form of the original 
Hamiltonian (j2ip . The Hamiltonian /i//,, preserves infinitely-many spaces 

pa,2,3) ^ (rf Va^Vf I < ni + 2n2 + 3^3 < n), n G N, 

with characteristic vector is (1,2,3), they form an infinite flag. The spectrum of hn^ is given 
by 

epi,P2,P3 = 2w(pi + 3p2 + 5p3), Pi = 0, 1, 2, ... , 
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with degeneracy pi + 3p2 + 5p3 = integer. It corresponds to the anisotropic harmonic oscillator 
with frequency ratios 1:3:5. Eigenfunctions cl)n,i of hn^ are elements of Vn''^'^\ The number 

(1 2 3) (12 3) 

of eigenfunctions in Vn ' is maximal possible - it is equal to dimension of Vn ' ■ 

(1 2 3) 

The space Vn ' is finite-dimensional representation space of a Lie algebra of differential 
operators which we call the h^^^ algebra. It is infinite-dimensional but finitely generated algebra 
of differential operators with 30 gener ating elements of 1'* (14), 2°'^ (10) and 3'^'^ (5) orders, 
respectively, plus one of zeroth order. They span 5-1-5 Abelian (conjugated) subalgebras of 
lowering and raising generators L and and one of the Cartan type algebra B (for details 
see [H]). The algebra B = gi2 © I2, where I2 is two-dimensional ideal. The Hamiltonian hn^ 
can be rewritten in terms of the generators of the /i^^^-algebra. 

By adding to hn^ (f2T]) the operator ([6]) in the variable ri (of the same type as for the Calogero, 
BCn and G2 models) we get the operator /i/^g +5h^'^°^^ which has the finite-dimensional invariant 
subspace 

Vk = {nP\0<p<k), 

of the dimension (A; -|- 1). Hence, this operator is quasi-exactly-solvable. Making the gauge 
rotation of this operator and changing variables r back to Cartesian ones we arrive at the quasi- 
exactly-solvable Hamiltonian of a similar type as for the Calogero, BCn and G2 models 



By adding to /ijjg ([2T|) the operator we preserve the property of exact-solvability. This 

operator preserves the flag and the emerging Hamiltonian has a form of (|2ip with the 

extra term ^ in the potential. Its ground state eigenfunction is (r^)'>'^o- It is the exactly- 
solvable generalization of the i^a-rational model (j2ip with the Coxeter group H3 as the discrete 
symmetry group, 

r 

^Ha = T-l-Hs + ^• 

2.10 Case H4 

The H/j^ rational Hamiltonian reads 



^ 5^ n , 1/(Z/ — 1) 



---2 ^u. -r -2" 



-|- 2l'(h' — 1) ; r ; 7 7 r- 

Jxi+ -iy^23;2+ -1^*3X3 + (-l)/^4a;4 2 

M2, 3, 4=0,1 

-h2i/(i/-i) y y - — . — ^ 1 — ^ ^, (22) 

where {i,j,k,l} = {1,2,3,4} and all even permutations, lo, v are parameters and 
¥'± = ^— , 

the golden ratio and it algebraic conjugate. Symmetry of the Hamiltonian (j22p is the H4 Coxeter 
group (the symmetry group of the 600-ce//). It has order 14400. In total, the Hamiltonian 
is symmetric with respect to the transformation 

Xi i — Xj , i — . 



^These subalgebras can be divided into pairs. In every pair the elements of different subalgebras are related 
via a certain conjugation. 
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The ground state function and its eigenvalue are 



A?^A^A^ exp ' ^0 = + 30u), 



where 



4 



Ai = i[xk, A2= n [2;i + (-ir^2 + (-irx3 + (-irx4], 

k=l A'2,3,4=0,l 

n n h+(-irv+2;i+(-irv-^fc+o-x;]. 

{i,j,k,l} A»i, 2=0,1 

Making a gauge rotation of the Hamiltonian 

and introducing new variables ri^2,3,4) which are invariant with respect to the H4 Coxeter group, 
in a form of polynomials in x of degrees 2, 12, 20, 30 (degrees of H4) , we arrive at the Hamiltonian 
in the algebraic form [15] 



where 



All = 4ri, ^12 = 24r2, A13 = 4OT3, Am = 6OT4, 
A22 = SSriTs + Srf T2, A23 = -4rfr| + lArfr^ - 8T4, 

38 8 
A24 = lOrfrl + 60riV2r3 + 40rf T4 - 600x3^, A33 = -ynr| + 28rf r2r3 - -riV4, 

A34 = 21QtItIt2. + 60riV2r4 - 180ri^ t| + 30t|, 

A44 = -2175riT|T3 - 450rf T|r4 - 1350tiV2t| - 600ti^T3T4, 

Bi = 8(1 + 30z^) - 4wTi, B2 = 12(1 + lOz^) rf - 2Aujt2, 

B3 = 20(1 + 6iy)TfT2 - 40wT3, B4 = 15(1 - 30z>)rfr| - 450(1 + 2z/)riV3 - 60a;r4, 

which is amazingly simple comparing the very complicated and lengthy form of the original 
Hamiltonian ()22p . It is easy to check that the algebraic operator hn^ preserves infinitely-many 
finite-dimensional invariant subspaces 

p(i,5,8,i2) ^ ^^ni r-2^n^^n, | q < „^ + + 8n3 + 12n4 < n) , n G N, 

all of them with the same characteristic vector (1,5,8,12), they form the infinite flag. The 
spectrum of the Hamiltonian hn^ (|23p has a form 

ekiMMM = 2uj{ki + 6/c2 + 10A;3 + 15/04), h = 0, 1, 2, ... , 

with degeneracy /ci + 6/^2 + 10A;3 + 15A;4 = integer. It corresponds to the anisotropic harmonic os- 
cillator with frequency ratios 1 : 6 : 10 : 15. Eigenfunctions (j)n,i of hfj^ are elements of 'p^'^'^'^'^\ 

fl 5 8 12) fl 5 8 12) 

The number of eigenfunctions in Vn ' ' ' is equal to dimension of Vk ' ' ' . 

The space T^n ' ' ' ^ is a finite-dimensional representation space of a Lie algebra of differential 
operators which we call the h^^^ algebra. It is infinite-dimensional but finitely generated algebra 
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B 



L- 



Pp{B) 



U 



Figure 3. Triangular diagram relating the subalgebras L, U and B. Pp{B) is a polynomial of the pth 
degree in B generators. It is a generalization of the Gauss decomposition for semi-simple algebras where 
p=l. 



of differential operators, with 183 generating elements of I''* (54), 2'^'^ (24), 3'''^ (18), 4''<^ (18), 
5""^ (28), 6""^ (5), 8""^ (20), 9""^ (1), 12''<^ (14) orders, respectively, plus one of zeroth order. They 
span 26 + 26 Abelian (conjugated) subalgebras of lowering and raising generators L and L'H and 
one of the Cartan type algebra B (for details see [l6] ) . The algebra B = gi2 © -^3 , where is 
three-dimensional ideal. The Hamiltonian hn^ can be rewritten in terms of the generators of 
the /i^^^-algebra. 

By adding to hn^ (p2|) the operator ([6]) in the variable ri (of the same type as for the Calogero, 
BCn and G2, models) we get the operator hn^ + 6h^'^°^\ which has the finite-dimensional 
invariant subspace 

Vk = {riP\0<p<k), 

of dimension (A; + 1). Hence, this operator is quasi-exactly-solvable. Making the gauge rotation 
of this operator and changing variables r back to Cartesian ones we arrive at the quasi-exactly- 
solvable Hamiltonian of a similar type as for the Calogero, BCn and G2 models jl5j . By adding 
to (p2]) the operator 47^ we preserve the property of exact-solvability. This operator 
preserves the flag ^(^5,8,12) ^i^g 

emerging Hamiltonian has a form of (j22p with the extra 
term ^ in the potential. Its ground state eigenfunction is (r^)'^^'o. It is the exactly-solvable 
generalization of the ff4-rational model (I22p with the Coxeter group H4 as the discrete symmetry 
group, 

r 

= 'Hh4 + -2- 



3 Conclusions 

• For rational Hamiltonians for all classical A^, BCn and exceptional root spaces G2, -F4, 
-£'6,7,8 (also trigonometric) and non-crystallographic -^3,4, l2{k) there exists an algebraic 
form after gauging away the ground state eigenfunction, and changing variables to sym- 
metric (invariant) variables. Their eigenfunctions are polynomials in these variables. They 
are orthogonal with the squared ground state eigenfunction as the weight factor. 

• Their hidden algebras are gl{N + 1) for the case of classical An, BCn and new infinite- 
dimensional but finite-generated algebras of differential operators for all other cases. All 
these algebras have finite-dimensional invariant subspace(s) in polynomials. 

• The generating elements of any such hidden algebra can be grouped into an even number 
of (conjugated) Abelian algebras Lj, Ui and one Lie algebra B. They obey a (generali- 
zed) Gauss decomposition (see Fig. [3|). A description of all these algebras will be given 
elsewhere. 

There are two solvable potentials in ID in [0,cx3) which can be generalized to D. 

^These subalgebras can be divided into pairs. In every pair the elements of different subalgebras are related 
via a certain conjugation. 
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3.1 General view (exact-solvability) 

ES-case generalization: 

2 2,7 , 22, 7(^) 

(a generalization by replacing the symmetry 0{N) by its discrete subgroup of symmetry given 
by Weyl (Coxeter) group). 

3.2 General view (quasi-exact-solvability) 

QES-case generalization: 

22,7, 6,1,4 , -2 2, 7(^) , 6 , I, 4 
UJ r H — K + ar + or — > u) r -\ 5 h ar + or 

(a generalization by replacing the symmetry 0{N) by its discrete subgroup of symmetry given 
by Weyl (Coxeter) group). 

In conclusion I must emphasize that the algebraic nature of all above-considered systems was 
revealed when: 

The invariants of the discrete group of symmetry of a system 
are taken as variables (space of orbits). 
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